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Abstract 

A mixing chamber used in rocket engine testing at the NASA Stennis Space Center 
is modelled by a system of two nonlinear ordinary differential equations. The mixer 
is used to condition the thermodynamic properties of cryogenic liquid propellant by 
controlled injection of the same substance in the gaseous phase. The three inputs of the 
mixer are the positions of the valves regulating the liquid and gas flows at the inlets, 
and the position of the exit valve regulating the flow of conditioned propellant. Mixer 
operation during a test requires the regulation of its internal pressure, exit mass flow, 
and exit temperature. A mathematical model is developed to facilitate subsequent con- 
troller designs. The model must be simple enough to lend itself to subsequent feedback 
controller design, yet its accuracy must be tested against real data. For this reason, the 
model includes function calls to thermodynamic property data. Some structural prop- 
erties of the resulting model that pertain to controller design, such as uniqueness of the 
equilibrium point, feedback linearizability and local stability are shown to hold under 
conditions having direct physical interpretation. The existence of fixed valve positions 
that attain a desired operating condition is also shown. Validation of the model against 
real data is likewise provided. 


1 Introduction 

The NASA John C. Stennis Space Center (SSC) conducts extensive ground-based testing 
and flight certification of rocket engines, in particular, of the Space Shuttle Main Engine 
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Figure 1: The Mixer Subsystem 

(SSME). Combustion chambers and turbomachinery related to rocket engines are also tested 
at SSC. This work is part of an on-going effort to develop a software package providing flex- 
ibility in simulation and control tasks [8, 2, 3] frequently found in test operations at SSC. 
Test conditions require that liquid propellants, namely liquid oxygen and liquid hydrogen 
(LH2) be supplied to the engine or component at very precise conditions of temperature, 
pressure and mass flow rate. An excess or deficiency in any of these three flow parameters 
may result in damaged components or in a sub-optimal test. To achieve the required con- 
ditions, the delivery system includes a mixing chamber, henceforth referred to as “mixer” . 
The mixer subsystem is depicted in Fig. 1. LH2 is stored in the run tank, which is kept 
at a constant pressure by an independent control loop. Gaseous hydrogen (GH2) is stored 
in high pressure bottles, and has a higher temperature than the LH2. One control valve, 
referred to as the “liquid valve” is used to manipulate the flow of LH2 into the mixer. 
There is also a “gas valve” and an “exit valve”. -The positions of the valves constitute 
the only control variables for the mixer subsystem. A mixer control system must achieve 
tracking and regulation of mixer outputs to desired values. The outputs of interest are the 
mixer pressure, exit flow temperature, and exit mass flow. A dynamic model and controller 
are intended to replace the current method of operation of the mixer, which uses only a 
steady-state thermodynamic model and extensive heuristics. The use of the modeling and 
control techniques described herein is expected to provide great flexibility and better mixer 
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performance. In this article, focus is placed on the mathematical modelling stage, together 
with an analysis of model properties and data validation. 

2 Transient Thermodynamics Model of the Mixer 

2.1 Valve Models 

The flow wi of LH2 through the liquid valve depends on the pressure difference across the 
valve, the density of the source LH2 and the liquid valve opening coefficient C v i according 
to [1]: 

= C\C vl V(Pl ~ P)Pi (1) 

In English units, Pi is the source LH2 pressure and P is the mixer pressure in psia; pi is 
the source LH2 density in lbm/ft 3 , and C\ = 1.76 x 10“ 2 . For gas flow, there is a separate 
expression [1]: 

C 2 C vg ^jPj - P 2 ,when P g < 2 P 

( 2 ) 

CiC vg -^ / Tgp ,when P g > 2 P (choked flow) 

where T g is the source GH2 temperature in degrees Rankine (°R), P g is the source GH2 
pressure in psia, and C vg is the gas valve opening coefficient. In English units, C 2 = 
2.857 x 10 -2 and C4 = 2.423 x 10~ 2 . Specifications for the conditioned exit flow require 
that the propellant be in the liquid phase. Therefore it is assumed that the form of Eq. (1) 
is also valid for the exit flow, namely 

w f = C 3 C vl y/(P - Ps)P ( 3 ) 

where p and P are the mixer density and pressure, respectively, and P s is the pressure at 
the outlet of the exit valve. 

2.2 Mass and Energy Balances 

The set of two nonlinear differential equations which constitute the model can be derived 
from the conservation of mass for compressible flow and from the First Law of Thermo- 
dynamics for a control volume including transient terms [9]. Upon manipulations, the 
equations may be expressed in terms of internal energy and density derivatives as: 

P = 
u — 


— (Wl + Wg - W e ) 

^ 7 1 u) -f- Wgi^Jig u) w e (h e r^)J 
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( 4 ) 

( 5 ) . 


where V is the fixed mixer volume, hi, h g and h e are the source liquid, source gas, and exit 
enthalpies, respectively. The mass flows Wj are functions of pressure, and are calculated 
according to Eqs. (1) and (2). 

2.3 State Variable Model 

Further manipulations are required to put the above physical model in state-space form. It is 
assumed that enthalpy is conserved across the exit valve, and the definition [9] h = u+CP / p, 
where C is a constant depending on the choice of units, is used. For English units, C = 0.185. 
The resulting model is compactly expressed as 

P = fi(P(p,u),p)C v i+ f 2 (P(p,u),p)C vg + f 3 (P(p,u),p)C ve (6) 

u = gi(P(p,u),p,u)C v i + g 2 {P{p,u),p, u)C vg + g 3 (P(p,u), p,u)C ve (7) 

where 


h = 
h = 
h = 
h — 

9i = 
g2 = 

gz = 


CW(P l -P) Pl 
C' 2 , if Pg > 2 P 

G'l'\j Pg ” P 2 ’ if^<2P 
cW(p-p*)p 


h 

h 

c'i 



The constants are given by C[ = C\/V, C' 2 — C 2 ^jTgp g /V , C' z = —CzjV = —C\/V, 
C'{ = CC' 3 and C[ = C^^T^Pg/V. It is to be noted that, given initial conditions of internal 
energy and density, the numerical solution of Eqs. (6, 7) requires that the mixer pressure 
P be known at all times as a function of p and u. Such computation poses difficulties 
and requires special routines. Indeed, strictly speaking, density and internal energy do not 
completely determine pressure. Thermodynamic data for a variety of substances shows that 
one may find density-energy-pressure triples that have the same densities and energies with 
distinct corresponding pressures. Although it is possible, in principle, to find correlations 
that link the variables in limited ranges, the number of separate expressions and their 
mathematical form makes this method inadmissible in the development of a controls model. 
Fortunately, the errors introduced by assuming that pressure is a function of density and 
internal energy are fairly small in the expected range of mixer operation. The selection 
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of a pressure based on p and u from thermodynamic data is not trivial, if correlations are 
not to be used. The available data consists of two tables, one giving density and the other 
energy, when pressure and temperature are known. The sought routines must perform then 
a reverse look-up of the tables. The details of how these routines work are out of scope, 
and it suffices to say that they have . been proved to be accurate by using them to obtain 
pressure and temperature from density and energy and then recovering input data by using 
the original tables, with little error. For further details, readers are referred to [5]. 

2.4 Output Definitions 

Mixer operation requires the simultaneous tracking of exit temperature, exit flow, and mixer 
pressure. Exit temperature and mixer pressure are functions of the energy and density 
states. Computation of these functions requires the intervention of thermodynamic tables 
and interpolation algorithms which cannot be represented in closed form. However [5], to 
each exit temperature and mixer pressure combination in the expected range of operation 
there corresponds a unique value of the state [p, u] T . Therefore, it is convenient to specify 
density and energy as outputs, along with exit mass flow. A real-time control system should 
perform pre-processing of commanded mixer outputs to obtain the corresponding desired 
state values. With these output definitions, the controls model becomes 


f 





' C* 
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Mp,u) h (p,u) h{p,u) 



< 
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9i(p,u) 92(p,u) 93 (p,u) _ 
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Cve 

„ j 


(</ = [a «, -Vf 3 c ve f 


where [C v i C vg C ve ] T is the control vector and y is the output vector. 

3 Model Properties 

3.1 Uniqueness of Equilibrium 

The equilibrium point indicates the steady values of the density and internal energy of 
the mixer when the control valves are set at fixed positions, for a given set of input flow 
properties. For given values of the input fluid properties and C v coefficients, setting p cv = 0 
establishes that m e = m g + rhi and results in an expression relating the equilibrium density 
p cv to the equilibrium mixer pressure P: 

Pcv = p(P) . ( 9 ) 
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Graphical Interpretation of Mixer Equilibrium 
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Figure 2: Mixer equilibrium 


Setting u cv = 0 gives the equilibrium exit enthalpy h e = h cv in terms of the input enthalpies 
and the mass flows: 


h cv = ^ + = h{ P ) 


( 10 ) 


m +m g 

The above enthalpy must match the thermodynamic property data at P and p cv of Eq. (9) 
that is 


h cv — hth{pcv,P) 


Substituting Eqs. (9) and (10) into the above equation results in a single expression which 
gives the equilibrium pressure: 

h(P)=h th (p(P),P) ( 11 ) 


A graphical interpretation of the equilibrium solution is shown in Fig. 2. The curve pP vs. 
P has been drawn on the base plane. This plane curve is mapped to a space curve by the 
thermodynamic property function h t h- The equihbrium exit enthalpy h(P) is a function of 
pressure only and the corresponding surface has also been graphed. The point where the 
space curve pierces the surface is the equihbrium point of the system. The monotonicity of 
the curve and surface indicate that there exists only one equilibrium point in the range of 
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interest of actual mixer operation. An iterative procedure is developed in [5] that calculates 
the equilibrium point, given thermodynamic parameters and fixed control inputs. 


3.2 Valve Positions for a Prescribed Operating Condition 


Having three independent controls in a two-state model allows the selection of steady values 
for the states, and, in addition, an extra degree of freedom is available. This degree of 
freedom can be used to fix the exit mass flow with the desired thermodynamic properties, 
as shown next. Supppose it is desired to have a given exit mass flow, with prescribed 
temperature (measured at the outlet of the exit valve). Let these quantities be denoted by w e 
and T e . The back pressure P s at the outlet valve exit is assumed to be constant. The problem 
is to determine the valve coefficients that achieve this. Two degrees of freedom are used for 
w e and T e , and the third one is used to meet a desired mixer operating pressure, P. The back 
pressure P s and T e determine the enthalpy h, which is constant across the valve. Therefore, 
the enthalpy and mixer pressure at the exit valve inlet are known and, in turn, determine 
the mixer density from thermodynamic data. A static energy and material balance gives 
the required input flows that achieve a prescribed flow-exit enthalpy combination. The 
following formulas are straightforward to derive: 


wi 


w 
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w e (h — h g ) 

hi hg 

WeijH ~ h) 

hi hg 


( 12 ) 


where h g and hi are the gas and liquid supply enthalpies, respectively. Using the exit flow 
and the input equilibrium flows from Eq.( 12), the three valve coefficients are determined 
from Eq. (1) and (2). This procedure is complementary to the one described in the previous 
section, which calculates the equilibrium point. 


4 Zero Dynamics and Feedback Linearizability 

Of the control techniques that are suitable for the present model structure, feedback lin- 
earization [4, 7] is the simplest and most direct. Input/Output linearization consists in 
finding a coordinate transformation which results in a system which is linear between the 
new inputs and the outputs. Linear controller design can then be applied to the transformed 
system. The technique, however, may have disadvantages. Specifically, the complexity in- 
troduced by the cancellation of all nonlinear dynamics may be to high for a realistic imple- 
mentation with limited computational resources. The other significant disadvantage is that 
all state measurements are required. Peaking of control signals is also a factor of concern, 
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especially when actuator saturation and rate limitations are present. Linearizability anal- 
ysis, however, is worth examining, since essential features of the model appear during the 
analysis and provide a basis for other controller design techniques. 

4.1 The Concept of Zero Dynamics 

A fair mathematical discussion of the concept of zero dynamics is out of scope and can be 
found in standard sources such as [4, 7]. Intuitively, when the order of the linear system 
arising from the input transformation is less than the order of the original system, the 
remaining zero dynamics can be realized by an appropriate set of state variables. These 
variables represent some physical aspect of the system and must be kept bounded. Unsta- 
ble zero dynamics constitute a fundamental limitation in several control design schemes, 
including feedback linearization. When a system has stable zero dynamics, it is said to 
be minimum-phase , a term which generalizes the concept of right-half plane zeros in linear 
systems. 

4.2 Relative Degree and Input Integration 

The partial relative degree [7] of an output of a system is the order of the lowest derivative 
of the output which is directly affected by at least one of the inputs. The model of Eqs. (8) 
has a characteristic which prevents direct application of input/output linearization theory. 

It is seen that the third component of y is algebraically related to one of the control inputs, 
namely C ve . This implies that the partial relative degree of y$ is zero. Conceivably, one 
could use a time-explicit control of the form 

n u v _ yu{t) 

C “ ( ’ A “>~ Vh[P(p,u),p) 

to attain perfect tracking of the output mass flow rate and use the two remaining controls 
to force the flow to have the desired pressure and temperature. This approach suffers from 
the drawback of not employing feedback and therefore of being not self-correcting or robust. 
One way to get around this problem is to augment the exit valve channel with an integrator, 
that is, let: 

where v is a new control input. Now C ve is regarded as a state, and the resulting system is 
of third order, with three inputs and outputs. If the arguments of functions Qi and fi are 
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dropped from the notation, the new system equations become: 

P = flCvl + }?p V g + fzC ve (13) 

u — 9iC v i + g 2 C vg + gzC ve (14) 

Cve = v (15) 

y = [p u -Vf 3 C ve f (16) 

Upon differentiating the outputs once, it is seen that the partial relative degrees are all 1: 

Vi — fiC vi + f 2 C vg + f 3 C ve (17) 

V 2 = Q\C v i + g 2 C V g + g 3 C ve (18) 

Vs = -v]vh + C ve (^^(f 1 C vl +f 2 C vg +fiC ve )+ 
d f dP 

+ (9i C v i + 9 c 2 .C vs + g 3 C ve ) 

+ §^C flCvl + hC vg + hc w ) ) ] (19) 

Upon rearranging, the output derivatives can be expressed compactly as 


y = D + Ew 

where y = [yi y 2 y 3 ] T , w = [C vi C vg v] T and 

fsC've 

D — I 9sGve 

~ VC ve[^(%f3 + ^9z) + ^afh] 

f 1 


E = 


9 1 


h 

92 


-VCvei^i^fl + ^9l) + §jy/l] -VC ve [^(^f 2 + §£ 52 ) + |^/2] 
Provided E is invertible in a region Q of the state space, the feedback law 

w = E- 1 (y d -F(y-y d )-D ) 


0 

0 

-Vfz _ 
( 20 ) 


achieves exact linearization of the system, with tracking error dynamics given by 

( V ~ Vd) + r (y -yd) = 0 


If r is chosen as a diagonal positive-definite matrix, the resulting control law is called “de- 
coupling control” , since the dynamics of the output errors are decoupled. If the appropriate 
function definitions are substituted, the forms of the E and D matrices are as follows: 


D = 


vc^cl 

2 


C3 \f (-P WPC W 
[p% + 0.185ff + (P-P.) 
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0 

0 


cW(Pi-P)pi h 

E= | ciV(Pi-P )fn (**=* ) h (^) 

— \/(A— -P)pi p VCW/2C3 p —VC / Tp W\fi 

■2^/{P-P,)p 1 -2^J (P—Pj ) ~p VC 3Vl.P P s)P 


where 


T; = 


r fl = 


' &P ,dP , 

w ,, ,ap 

- ^)~T~ + (-P - ^s) 


<9p 


4.3 Input-Output Linearizability of Augmented Model f „ 

The ability to construct a feedback linearization controller hinges, first, on stable zero 
dynamics of the augmented system, and, second, on the inyertibility of matrix E. 


4.3.1 Invertibility of E 


By inspection, it is readily seen that the first two rows of E are linearly independent provided 
hg ^ hi. This has a direct physical interpretation: if the two fluids have the same thermal 
properties (he., enthalpies), the ability to change the thermal properties of the mixture by 
changing the relative flows is lost. Gas and liquid enthalpies are different for the expected 
operating conditions. The third row is linearly independent from the first two if P — P s > 0, 
which is also true for the mixer. Therefore E is invertible ovef the whole range of expected 
mixer operating conditions'. * * * ^ ■ > 5 , t ‘ 

'■ » x i 

4.3.2 Zero Dynamics of the Augmented Model 


As it is known, the zero dynamics is preserved under input transformations [4]. This implies 


that one may examine the non-augmented model for zero dynamics and draw conclusions 
about the augmented model’s zero dynamics. Suppose the exit flow is to be held constant 
at a value Tjo- The only way in which this can be achieved is by letting 


P've — 


*30 

■Vh(t) 


( 21 ) 


at all times. Differentiating the other two outputs and equating them to zero results in * 


C v i{t) — + -Ej (22) 

°’ ,(t) = 'ik( siWC '" w+ ^fi 5 ) (23) 

Upon substitution and rearrangement, it is seen that there exist three uniquely defined 
control inputs which hold the outputs constant. Since two of the outputs coincide with 
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the states, it trivially follows that they are kept bounded, and therefore the system has 
stable zero dynamics (i.e., the system is minimum-phase). The control inputs are given by 
Eq. (21) and, dropping the time notation, 


Ci 


vg 


C v i 


^30 A f 33 _ 2i 

v U h) 


fl 92 - 91/2 


y 30 


1 - 


/ 1/2 


Vfi L fl92 - 91 f 2 




(24) 

(25) 


Note that the above formulas can be used to find the valve positions at which the system 
has a prescribed outflow and a pair of thermodynamic properties. 


5 Small-Signal Model 

The dynamic model of the mixer can be written in the form 

%i — F\(zi, Z 21 C V g , C v i , C ve ) (26) 

*2 = Fz(z 1 , z 2i C vg , C vh C ve ) (27) 

where fq (•) and 7*2 (•) are nonlinear functions of the state and valve coefficients. For the 
remainder of this article, we will denote constant or equilibrium values of any variable by 
an upper bar (•). Given constant values of valve flow coefficients C v = [C vg C v i C ve ] T 
(superscript T denotes transposition) and constant fluid properties, the state of the model 

z\ : Internal Energy 

Z 2 : Density 

reaches a constant equilibrium point z. Next, consider perturbing such an equilibrium by 
small signals x(t) and u(t) so that 

z(t)=z + x(t) and C v (t) = C v + u(t) 



where u(t) denotes a small, valve-coefficient correction/regulation signal. Then, a standard 
linearization of Eqs. (26)-(27) results in the small signal model 

x = Ax + Bu (28) 

where x(t) is the small perturbation state vector, u(t) is the small perturbation control 
signal, and the two-by-two matrix A and two-by-three matrix B are given by 



" dFi 

dF± 


‘ dFj 

dF± 

dFi 1 


dzi 

dz2 

; 5=j 

dC v g 

dc vl 

dC ve 


dF 2 

dF 2 

1 dF 2 

0F 2 

dFz. 


_ dzi 

dz 2 j 


. OCvg 

dC v i 

dC ve . 
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where the partial derivatives are evaluated at the equilibrium state z and constant valve 
flow coefficient vector C v . Output equations of the form 


y = Cx + Du 


where matrices C and D are appropriately dimensioned are easily appended to the model (28) 
to account for the measurement of certain variables such as temperature, pressure, or flow. 
An output of interest is the exit flow w e . Using Eq. (3) and linearizing around a chosen 
equilibrium point, we obtain the linear approximation 

da e da e 


W elin 


dzi dz2 


x + [a e \ 


eq 


eq 


U v 


— C'vflow'E “f” 3D flow'U've 


Other outputs of interest axe the mixer pressure P v and the exit temperature T e . Taking 
these to be functions of the mixer internal states, that is, 


Pv = and T e = f T (zi(t),z 2 (t)) 


then, the linearization around the chosen equilibrium gives 

~dP v dP v 


Pylin — 


dzi dz 2 


x — CpyX 


eq 


and 


Pelin — 


X = C te X 


~ dT e 8T e 
dzi dz 2 

The indicated partial derivative terms may be calculated numerically using accesory rou- 
tines. 


eq 


5.1 Numerical Example 

Table 1 lists the data for a typical equilibrium point EQi corresponding to a given set of 
outflow requirements. The matrices A and B are evaluated at the equilibrium EQi listed 
in Table 1 to be 


-17.52 

-3393.08 " 

; b = 

’ 25.11 

332.74 

-42.01 

-0.1597 

-43.45 


0.651 

0.6612 

-0.4919 

£ - 


For an output y = w e , the linear approximation is 

Welin = C v f low x + Dfi ow u ve = [0.2386 66.84] x + [1.23 ]u ve 

The eigenvalues of A he in the left half of the complex plane, establishing the local stability 
of the model by Lyapunov’s linearization method [10]. Although not shown here, the 
computation can be performed at any point in the range of operation, showing stability in 
a larger region. Local controllability [10] can be likewise evaluated from the linearization at 
various points in the operating region. In this example, the controllability matrix [. B | AB] 
has rank two. 


12 



1 

P 

(psi) 

T 

(F) 

P 

{lbm/ft 3 ) 

h 

{BTU/lbm) 

Cy open 

%OPEN 

w 

Ibm/s 

Zi 

*2 

gh 2 in 

m 

90 

2.91 

2113.6 

230 

0.96 

- 

- 

- 

gh 2 out 

6804 


' 1.702 

2113.6 .. 


- 

3.72 

- 

** 

lh 2 in 


-340 

! 5.042 

329 

115 


- i 

- 

- 

lh 2 out 


-326 

4.519 

329 


1 

33.21 


- 

Mixer 

6804 

-273.32 

3.84 

508 j 


- 

- 

181.39 

3.84 

Outflow 

5533 

-266.21 

3.404 

508 1 

270 

11.12 

36.93 "| 

- 

- 


Table 1: Initial Equilibrium Data. Underlined Values are Computed. (IN=Into Valve; 
OUT=Out of Valve). z\ in BTU/lbm; z 2 in lbm/ft 3 . 

6 Model Validation 

The usefulness of the model must be demonstrated with real data. Visual and quantitative 
assesments of the model’s ability to capture the essential dynamics of mixer operation axe 
provided in this section. Data available for model validation consisted of a preliminary 
operation of the mixer with Nitrogen as working substance. The thermodynamic lookup 
routines discussed earlier were set to work with Nitrogen, without modifications to the 
model equations. 

6.1 Calibration of Valve Coefficients 

Part of the experimental data corresponds to nearly steady conditions. A measurement 
of the flows and of some thermodynamic properties can then be used to estimate the cor- 
responding values of the model valve coefficients. The coefficients corresponding to other 
percent openings are found by direct proportionality. 

6.2 Boundary Conditions and Valve Position Histories 

The thermodynamic properties at the inlet of the liquid and gas valves change with time 
and were measured throught the duration of the experiment. In particular, pressure and 
temperature were recorded, data from which the enthalpy and density can be easily found 
by direct lookup of appropriate thermodynamic tables. The exit fluid was dumped into 
the atmosphere, thus establishing a fixed boundary pressure at the outlet of the exit valve. 
Knowledge of the above properties, valve position history and initial conditions of the mixer 
states is sufficient to obtain a numerical solution to the differential equations. The initial 
conditions are obtained from the measured mixer internal temperature and pressure at the 
begining of the experiment. These mixer variables were measured at all subsequent times 
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Figure 3: Inlet Thermodynamic Properties and Valve Position History 

during the experiment and provide the basis for comparison and evaluation of the model. It 
is to be noted that during the actual experiment, the exit valve was kept at a fixed position. 
Figure 3 shows the measured temperatures and pressures of the liquid and gas nitrogen as 
recorded, as well as the valve positions. 

6.3 Simulation 

The boundary variables shown in Figure 3 were fed into a Simulink model of the mixer 
equations. Figures 4 and 5 show the simulated and measured temperature and pressure 
inside the mixer. Despite of some spurious effects due to the interaction of the integration 
method with noisy input data, it can be seen that the agreement is very good. It should be 
noted that the pressure curve has been slightly shifted to compensate for an initial offset 
error found in the original data. 

6.4 Immediate Improvements to the Model 

The mathematical model developed may be improved and extended by considering the static 
and dynamic characteristics of valve actuation. In the above development, it is assumed 
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Figure 5: Pressure Validation 
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that the valve opening coefficients can be directly and instantaneously commanded to any 
desired value. This assumption did not interfere, however, with the model validation, since 
actual valve coefficient responses were available and were used instead of valve commands. 
In the practical situation, the control system has authority over the commanded openings 
only. A simple model of valve response is given by the first order system 

C v = --C v + /3( (29) 

T 

where ( is the commanded valve opening, r is a time constant and /3 is a proportionality 
constant. An equation of the above form must be written for each of the three valves in the 
system and the model must be augmented with these dynamics, thus becoming a 5th-order 
model. Note that direct proportionality has been assumed to exist between ( and C v under 
steady conditions ( C v = 0). If nonlinearity is observed in an experimental characterization 
of the valve, it can be included in the form: 

C v = --C v +g{0 (30) 

T 

where g is a function to be determined from static calibration, as done in [6]. Further 
improvements to the model include adding a disturbance term -mainly corresponding to 
unmodelled heat transfer- to be used in robust controller design. A Matlab-based graphical 
user interface (GUI) is currently under development that will provide flexibility in using the 
model in simulation studies. 

7 Conclusions 

A mathematical model for the mixer system is presented. The model consists of a system 
of two nonlinear equations having density and internal energy as states. The independent 
variables are the valve positions represented by flow coefficients, and the controlled outputs 
are mixer pressure, exit temperature and exit mass flow. Model reliability is ensured by the 
use of function calls to real thermodynamic data. It is shown that, in the expected range of 
operation, the mixer has a single equilibrium point for each set of fixed valve coefficients and 
that it is feedback linearizable and locally stable. The reachability of operating conditions 
is shown by the possibility of determining a unique set of valve coefficients that achieve 
a desired set ot steady outputs. The model is validated against experimental data with 
excellent agreement. 
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